Anomalous B meson mixing and baryogenesis 
in a two Higgs doublet model with top-charm flavor violation 
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There exist experimental hints from the B sector for CP violation beyond the Standard Model 
(SM) CKM paradigm. An anomalous dimuon asymmetry was reported by the DO collaboration, 
while tension exists between B — > rv and S^k- These measurements, disfavoring the SM at the 
~ 3a level, can be explained by new physics in both B®- B® and B°- B® mixing, arising from (1) new 
bosonic degrees of freedom at or near the electroweak scale, and (2) new, large CP-violating phases. 
These two new physics ingredients are precisely what is required for electroweak baryogenesis to work 
in an extension of the SM. We show that a simple two Higgs doublet model with top-charm flavor 
violation can explain the B anomalies and the baryon asymmetry of the Universe. Moreover, the 
presence of a large relative phase in the top-charm Yukawa coupling, favored by B® s -B® s mixing, 
weakens constraints from ejc and b — s- 57, allowing for a light charged Higgs mass of 0(100 GeV). 
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I. INTRODUCTION 

Precision tests of CP violation have shown a remark- 
able consistency with the Standard Model (SM), where 
all CP-violating observables are governed uniquely by the 
single phase in the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix PJ. Yet the search continues. Many well- 
motivated extensions of the SM, such as supersymmetry, 
contain new sources of CP violation at the electroweak 
scale. Furthermore, new CP violation beyond the CKM 
phase is likely required to explain the origin of the baryon 
asymmetry of the Universe. 

Recent analyses have suggested that the CKM 
paradigm may be in trouble. First, the DO collabo- 
ration has measured the like-sign dimuon asymmetry, 
arising from CP violation in the mixing and decays of 
B® mesons, in excess over SM prediction at the 3.2(7 
level 0. Second, there is tension at the ~ 3er level be- 
tween the branching ratio for B + — > t + v and the CP 
asymmetry S^k in B% —t J/ip K 0, Q • Additionally, 
CDF and DO have measured the CP asymmetry 
in B® — > J/ip(f>. While their earlier results (each with 
2.8 fb -1 data) showed a ~ 2a deviation from the SM 
this discrepancy has been reduced in their updated anal- 
yses with more data (5.2 and 6.1 fb -1 , respectively) @. 

Although further experimental study is required, taken 
at face value, these anomalies suggest CP violation from 
new physics (NP) in the mixing and/or decay amplitudes 
of B® and B® mesons Q. Recently, the CKMfittcr group 
has performed a global fit to all flavor observables, al- 
lowing for arbitrary new physics in B® S -Bj s mixing am- 
plitudes 0. They conclude that the SM is disfavored at 
3.4cr, while the data seem to favor NP with large CP- 
violating phases relative to the SM in both B® and B® 
mixing. At the level of effective theory, this NP takes 
the form 
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These operators can arise from new bosonic degrees of 
freedom at or near the weak scale, with new large CP- 
violating phases (914121] . 

It is suggestive that the same NP ingredients, new 
weak-scale bosons and new CP violation, can also lead 
to successful electroweak baryogenesis (EWBG). EWBG, 
in which the baryon asymmet ry is g enerated during the 
electroweak phase transition |13H15| . is particularly at- 
tractive since two out of three Sakharov conditions [l6| 
can be tested experimentally First, a departure from 
thermal equilibrium is provided by a strong first-order 
phase transition, proceeding by bubble nucleation. While 
this does not occur in the SM [l7| , additional weak-scale 
bosonic degrees of freedom can induce the required phase 
transition; these new bosons can be searched for at col- 
liders. Second, there must exist new CP violation be- 
yond the SM [18[. This CP violation must involve par- 
ticles with large couplings to the Higgs boson, since it 
is the interactions of those particles with the dynam- 
ical Higgs background field that leads to baryon pro- 
duction. Precision tests, such as electric dipole moment 
searches [l9[ and flavor observables, can probe directly 
CP violation relevant for EWBG. (The third condition, 
baryon number violation, is provided in the SM by weak 
sphalerons [2[j; however, its rate is highly suppressed in 
all processes of experimental relevance.) 

If we wish to connect Eq. ((T|) to EWBG, it is bet- 
ter to generate these operators at one-loop, rather than 
tree-level. Constraints on the mass differences AMd. s in 
the B% s systems require that A 2 /|c d | > (500 TeV) 2 and 
A V|c s f > (100 TcV) 2 (HI- For tree-level exchange, it 
seems unlikely that all three Sakharov conditions can be 
met at once. Sufficient baryon number generation typ- 
ically requires couplings > O(10 _1 ), such that Cd, s ^ 
0(1O~ 2 ), while a viable phase transition requires A < 1 
TeV. Therefore, EWBG requires A 2 /|c d , s | < (10 TeV) 2 , 
at odds with AMd, s constraints. However, if the opera- 
tors in Eq. ([T]) arise at one-loop order, Cd. s will have an 
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additional l/(4ir) 2 loop suppression, allowing for both 
large couplings and lighter scale A, without conflicting 
with AMd, s constraints. 

In this work, we propose that a simple two Higgs dou- 
blet model (2HDM) can account for both anomalous CP 
violation in a -&] mixing and EWBG. Previous works 
have studied CP violation in B® S -B® s mixing within a 
2HDM [Ul^. Our setup, described in SecHU is different: 
we assume the NP Higgs doublet (H + , H° + iA°) medi- 
ates top-charm flavor violation. In this case, the NP - 

B® mixing amplitudes (M^^np are generated at one- 
loop order through char ge c urrent interactions mediated 
by H + (similar to Ref. |12|). rather than through tree- 
level exchange [9j4TT|. In Sec. IIII1 we compute {M^^np 
in our model. We find: 

• The best fit values to both Mf 2 and Mf 2 , from 
Ref. Q, can be explained in terms of a single NP 
phase d tc (defined below). 

• For large values of i? tc prefered by B® S -B® s mixing 
observables, constraints from ck and b — > sj are 
weakened and can be light (m#± ~ 100 GeV). 

In Sec. [TVl we discuss in detail EWBG in our 2HDM 
model. We focus on the CP violation aspects of EWBG, 
computing the baryon asymmetry in terms of the un- 
derlying parameters of our model by solving a system 
of coupled Boltzmann equations. We find that the pa- 
rameter region favored by flavor observables (specifically, 
a large iptLH® coupling) can easily account for the ob- 
served baryon asymmetry. However, the relevant CP- 
violating phase is unrelated to the phase i?t c entering 
flavor observables. In Sec. [VJ we summarize our conclu- 
sions. 



II. MODEL 



where the left-handed quark doublet is Ql = {ul, VcLl)- 
The SU(2) L contraction is HiQ L = H+(Vd L ) - H°u L . 
The 3x3 Yukawa matrices yu and yu couple right- 
handed it-type quarks ur = (u,c,t)n to left-handed 
quarks uj, = (u, c, €)l and = (rf, s,6)i- Working in 
the mass eigenstate basis, the matrix 

yu = di&g(y u ,y c ,y t ) = diag(m u , m c , m t )/v (4) 

is a diagonal matrix of SM Yukawa couplings, and V is 
the CKM matrix. Analogous Yukawa couplings arise for 
down quarks and charged leptons: 

J^yuk D (5) 

- MvdHI + y D H\)Q L - e R {y L H\ + y L H\)L L + h.c. 

where y D = dmg(y d ,y s ,y b ) and y L = diag(y e , y^, y T ) are 
the SM Yukawa couplings. 

The NP Yukawa matrices yu,D,L can be arbitrary. 
However, the absence of anomalously large flavor- 
violating processes provides strong motivation for an or- 
ganizing principle. In this work, we assume that flavor 
violation arises predominantly in the top sector. Specifi- 
cally, we take 

/ \ 

Vu = _0 , y D , L = . (6) 
\ y tc yu ) 

That is, we consider a hierarchical structure where the 
tR-tL and tR-CL couplings are dominant (with \ytt\ 3> 
|2/t c |), while others are suppressed. The zeros in Eq. ([6]) 
are meant to indicate these subleading couplings that for 
simplicity we neglect in our analysis. In our setup, flavor 
violation in meson observables arises at one-loop order 
through H ± charge current interactions, discussed in the 
next section. 



In a general (type III) two Higgs doublet model |22| . 
where both Higgs fields couple to each SM fermion, one 
can perform a field redefinition such that only one Higgs 
field acquires a real, positive vacuum expectation value 
(vev) [23|]. We denote the two Higgs doublets by 



III. FLAVOR CONSTRAINTS 



G+ \ 




( H+ \ 


h°+iG° 1 


, H 2 =\ 


H°+iA° 1 


V2 J 




V V2 J 



(2) 



where h°, H° (A ) are the neutral (pseudo)scalars, iJ is 
a charged scalar, and G ±,cl are the Goldstone modes eaten 
by the electroweak gauge bosons. The vev is v w 174 
GeV. In general, the physical neutral states can be ad- 
mixtures of h°, H°, A , depending on the details of Higgs 
potential. We neglect mixing in our analysis; in this case, 
Hi is exactly a SM Higgs doublet. 

The most general Yukawa interaction for u-typc quarks 

is 



^yuk 3 u R {yi/Hi + yuH 2 )QL + h.c. 



(3) 



Mixing and CP violation in the B®-B® system {q = d, s) 
is governed by the off-diagonal matrix element M\ 2 ~ 
|r? 2 in the Hamiltonian 0,IH, with Mf 2 (rf 2 ) asso- 
ciated with the (anti-)Hermitian part. Only the relative 
phase 4> q = arg(— Mf 2 /rf 2 ) is physical. The relevant ob- 
servables are the mass and width differences between the 
two eigenstates 

AM, = 2|M« 2 | , Ar g = 2|r? 2 | cos0 9 , (7) 

and the wrong sign scmilcptonic asymmetry 

, _ £Cgg -> fi+X) - Y{B° q -> yrX) _ \y\ 2 \ . 
° d - r(B° -> ^X) + T(B° -+ p-X) ~ \Mf 2 \ Sm<P " ■ 

(8) 

The dimuon asymmetry measured by DO arises from 
wrong sign semileptonic decays of both 59 and B® 
mesons and is given by A b sl 0.5 + 0.5 a s sl 
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d.L,s L - 



W ± ,G ± ,H ± 
—- — d L ,s L 



are 



FIG. 1: New physics B" d -B% and mixing amplitudes 

(M^ s )np arising from box graphs with exchange. 

In the SM, the mixing amplitude Mf 2 arises from 
box graphs, while the T\ 2 comes from tree-level decays. 
Therefore, it is plausible that NP effects enter predomi- 
nantly through mixing. Deviations in Mf 2 from the SM 
can be parametrized by 



Mf 2 = (Mf 2 ) SM + (Mf 2 ) NP 



(M? 2 ) SM A q . (9) 



The consistency of AMd, s with SM predictions constrains 
|A rfjS | w 1, at the 0(20%) level [|, while the dimuon 
asymmetry measurement disagrees with SM prediction 
at 3.2cr and requires 0(1) NP phases <j)f = arg(A 9 ) 0. 
Phases <fi^ also enter into CP asymmetries due to inter- 
ference between B® s decay amplitudes with and without 
mixing: e.g., the asymmetry for — > J ftp Kg is S^k s = 
sin(2/3 + <ft), with CKM angle /? = arg(-V cd V c * b V t * d V tb ) . 
As emphasized in Ref. Q, the presence of non-zero cj>^ 
can alleviate tension between k s and Br(_B + — > t + v) 1 
which is sensitive to j3 but not 4>^. 

To quantify these tensions, the CKMfitter group per- 
formed a global fit allowing for arbitrary A c i. s (dubbed 
"Scenario I"), findingthat the SM point (A d = A s = 1) 
is disfavored at 3.4cr Moreover, their best fit point fa- 
vors NP CP-violating phases in both B^ and B® mixing: 



(-12 



-3.3\ 
-3.4J 



and 



(-129 



-t^)°U(-51.6™) 

observables pre- 
dominantly through mixing, via box diagrams shown in 
Fig. □ We find 2 



51 _ 6 +i4.i^o 1 



In our model, NP effects enter B° 



q — 1 + C bq F 1 (xH,X t )/S (x t ) +C b F 2 {xH,Xt)/S {Xt) 



A 

where 



(yuV)u(SuV)* tj 

4V2G F m 2 w VuV* ' 



(10) 
(11) 



The tfid l L H + charge current couplings are (yuV)u = 
VttVu + VtcVci, f° r i = d,s,b. The NP loop functions 



1 Ref. Q did not include in their fit updated CDF and DO results 
for 0|, which showed improved consistency with the SM 
over previous results favoring non-zero <f>^ [EJ. 

2 We neglect running between the scales rat, Tnyvt an< f rn H , inte- 
grating out these degrees of freedom at a common electroweak 
scale. Moreover, we have neglected a NP QCD correction factor 
rj(x[i,xt)/riB arising at next-to-leading order [26t . 



F 1 (x H ,x t ) 



X t X H {xH - 4) logSg _ x t {x t - 4) 
(x H - l)(x H - %t) 2 (x t - 1)(xh — X t ) 
x t (x H x 2 - 2x H x t + Ax H - 3x 2 ) \ogx t 



F 2 (x H ,x t ) =■ 



(x t -l) 2 (x H ~x t ) 2 
c 2 - 2x t x H log(x H /x t ) 



(x H - Xt) 3 



(12) 
(13) 



where xt,H 



' H± /m 2 w , and S (x t ) « 2.35 is the SM 



loop function (e.g., see [25] ). 



-B\ s mixing from box graphs in a 2HDM have 
been computed previously [2?j . Here, a novel feature 
arises from the NP CP-violating phase associated with 
Vtc HI- We can write (yuV) u as 



(yuV)tb — yttVtb 



(yuV)td = yttVtd 1 



0+ 


ytcVcs 


yuVts 




ytcVcd 


yttVtd 



,vdtc 



(14) 



where tf tc = arg(y tc F cs ?/* t F t *). In the limit \y a \ > \y tc \, 
we neglect the term ytcVcb for i = b; however, y tc is non- 
negligible for i — d,s because the y u terms are Cabibbo 
suppressed. 

The NP phase that enters (M( 2 )np is i)tc, while for 
(Mf 2 ) NP it is (i?i C + /?), due to the different CKM struc- 
tures of (yuV)tB and (yuV)td- The best fit values for 
4>^ s are quite different numerically, but due to this extra 



, we can explain both 



phase "dtc 
(M! 2 ) 

as (Mf 2 ) S M.) 



(For y tc = 0, our model gives 



in terms of the single NP 

= 0, since 
2 



d . s 

p would have the same complex phase {VtbV t * q ) 



Our results for B® g - 



B° 



mixing are shown in Fig. [2] 
Here, we map best fit regions for Ad, s from Ref. [|[ into 
the parameter space of our model. We fix \y tt \ and m H ± 
and evaluate the prefered regions for |yj c | and "dtc con- 
sistent with B® S -E>d s mixing constraints. (As discussed 
below, EWBG favors \y tt \ ~ 1 and m H ± < 500 GcV.) 
The blue (red) contours correspond to the best fit re- 
gions at la (inner) and 2a (outer), for A^ (A s ). Since 
Ad, s are quadratic functions of |y tc |e i1?tc , the best fit re- 
gions for Ad, s each map into two best fit regions in \yt c \, 
'dtc parameter space. 

We also implement constraints on our model from b — > 
s"/ and €k ■ The branching ratios for b — > S7, as measured 
experimentally [2^ and evaluated theoretically in the SM 
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FIG. 2: Top- charm flavor violation parameter space (\ytc\, #tc 
68% and 95% CL regions for Ad (A B ) from Ref. \Q] shown 
contours is consistent with ck at 68% (95%) CL. Light (dark 

at next-to-leading order (NLO) [30| . are given by 3 : 

BR[B X^]™^ 6GcV = (3.55 ±0.24 ±0.09) x 10~ 4 

(15) 

BR[B X, 7 ]gf >1 . flGeV = (3.60 ± 0.30) x 10~ 4 . 

We evaluate SM+NP contributions to BR[B -> A s7 ] in 
our model at NLO following Refs. [3(| HH, except that 
we take as inputs the best fit CKM parameters given 
in Table 11 of Ref. Q. Adding all errors in Eqs. (|T5|) 
in quadrature, we take the following constraint on our 
model: 

BR[B -> JT.7]|^GeV = ( 3 - 55 ± °- 39 ) x 10 ~ 4 ■ ( 16 ) 

In Fig. [21 the white (light grey) region corresponds to 
\ytc\i $tc parameter space consistent with Eq. (jTHJ) at less 
than lcr (2a), while the dark grey region is excluded at 
2a. 

NP contributions to K°-K° mixing arise in our model 
through box graphs analogous to Fig. [1] The strongest 
constraint is due to ex- In the SM, \ck |sm = (1-90 ± 
0.26) x 10" 3 [11, while experimentally \e K L P = (2.228± 
0.011) x 10-3g3l- The SM + NP value of £k is 

|eif|sM+NP = K e C e B K lm [{V ts V t * d ) 2 r]2 

X (Sofa) + C a dFx{xH,X t ) + c 2 sd F 2 {x H ,x t )) 

+(v cs v; d ) 2 m So(x c ) + 2(v cs v; d v ts v t * d ) m s (x c , Xt )] , 

(17) 

where NP enters through the coefficients c sd defined in 
Eq. dHJ). (We neglect NP NLO corrections to ?y 2 .) The 



3 In the observed value, the first error is experimental, while the 
second is a theoretical error associated with a photon shape func- 
tion used to extrapolate the branching ratio to different photon 
energies Ej. Also, although BR[_B — > X s ^\ has been computed 
at NNLO in the SM [H, wc work at NLO since 2HDM contri- 
butions have been computed at NLO only. 



consistent with flavor observables, for two choices of\y~tt\, rn,H±. 
by blue (red) contours. Region within dark (light) dashed green 
grey region is excluded at 68% (95%) CL from BR(B — > X 3 j). 

remaining SM input parameters in Eq. (|T7| are defined 
and tabulated in Ref. Q. Assuming a theoretical error 
bar as in Ref. [33j ] . we take the following constraint on 
our model 

Msm+np - (2.23 ± 0.30) x 10" 3 . (18) 

It appears that since |eif|sM < |e_R"|exp, this constraint 
would favor a small, positive contribution from NP. How- 
ever, \ck |sm itself is shifted to a central value \ck |sm = 
2.40 x 10~ 3 because the best fit CKM parameters in the 
presence of NP in B d S -B d s mixing (given in Table 11 of 
Ref. ||) are different than in a SM-only fit. As a result, 
Eq. (jT5J) favors a small, negative contribution from NP. 
In Fig. [21 the parameter region within the dashed dark 
(light) green contours is consistent with ck constraint in 
Eq. (TTHD at lcr ( 2cr )- 

Here, we make several important points. 

• Despite the fact that <fi d and <j>^ are quite differ- 
ent numerically, there exists regions of parameter 
space where both NP in B° d -B° d and B°-B° can be 
explained by a single phase $ tc . The la best fit re- 
gions for Ad >s overlap within the parameter space 
of our model (neglecting correlations between A^ 
and A s ). 

• The A s region that overlaps with the A^ region 
in Fig. [2] corresponds to the (f>£ = (— 51.6±9 4 4 1 )° 
solution. Therefore, our model predicts AT S > 0. 

• Although b — > s 7 and ck constrain a large para- 
metric region of our model, these two observables 
are consistent with observation in regions favored 
by B mixing observables. 

• A large phase § t c can weaken b — > s 7 and ex con- 
straints, and a light charged Higgs (m H ± ~ 100 
GeV) is not excluded. 

• The values of (\yu\y m H ± ) shown in Fig. [5] are con- 
sistent with R b = BR[Z -> bb]/BR[Z -> hadrons] 
at 95% CL [13 . 
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Although we chose only two illustrative values 
(\y u \,m H ±) = (0.8, 100 GeV) and (1.2, 350 GeV) in 
Fig. [21 there exists a consistency region between all these 
observables for parameters \yu\ ~ 1, \vtc\ ~ 0.05 — 0.1, 
and d tc ~ 3tt/4, for 100 < m H ± < 500 GeV. As we 
discuss below, EWBG favors \ytt\ ~ 1 and m H ± < 500 
GeV. 



IV. ELECTROWEAK BARYOGENESIS 

Given a NP model, viable EWBG requires: (1) the 
electroweak phase transition must be strongly first or- 
der to prevent washout of baryon number, and (2) CP 
violation must be sufficient to account for the observed 
baryon-to-entropy ratio Y£ bs rts 9 x 10 -11 . EWBG in 
a 2HDM has been studied many times previously [35| . 
Most recently, Ref. [36[ showed that a strong first or- 
der phase transition can occur in a type-II 2HDM for 
m h o < 200 GeV and 300 < m H o < 500 GeV. Although 
our 2HDM is not exactly the same as in Ref. J3f|, we 
assume that a strong first order transition does occur. 
(The phase transition can also be further strengthed or 
modified by the presence of scalar gauge singlets [37j or 
non-renormalizable operators [Hl|.) 

We now study baryon number generation during the 
phase transition. The dynamical Higgs fields during the 
transition gives rise to a spacetimc dependent mass ma- 
trix M(x) for, e.g., w-type quarks: 

i^mass = -UrMul + h.c. , M = y v vi(T) + y v v 2 {T) 

(19) 

where vi^iT) = (Hi 2 )t=£o are the vevs at finite tem- 
perature T ss 100 GeV. At zero temperature, when 
vi(T), V2{T) — > v, 0, we recover the usual T = masses. 
However, if i)2(T) 7^ 0, then CP-violating quark charge 
density can arise from yjj, as we show below. Left-handed 
quark charge, in turn, leads to baryon number production 
through weak sphalerons. In previous studies, CP asym- 
metries were generated by a spacetime-dependent Higgs 
vev phase, arising from CP violation in the Higgs sec- 
tor [3a, [36| . Here, we assume that the Higgs potential is 
CP-conserving, such that vi^{T) do not have spacetime- 
dependent phases and can be taken to be real. 

Is it plausible that V2 (T) ^ during the phase transi- 
tion? Following [l(J, the most general potential for Hi 2 
can be written 

V = A(F 1 t i/ 1 - v 2 ) 2 + m 2 H2 HlH 2 + XiH[H x hIH2 
+ X 2 h\H2HIHi + [X z {h\h 2 ) 2 + XiH\H 2 H\H2 
+ X h HlHi{H\Hi - v 2 ) + h.c] + X 6 {hIh 2 ) 2 (20) 

Our basis choice that (H 2 )jiii = requires that no terms 
linear in H2 survive when Hi -> v. The same state- 
ment does not hold at T ^ due to thermal correc- 
tions to V. First, since we expect i>i(T) 7^ V, terms 
linear in H2 appear proportional to A5. Second, top 



quark loops generate a contribution to the potential 
(ytyttT 2 HIH2 1 '4 + h.c), given here in the high T limit, 
also linear in H 2 - A proper treatment of this issue re- 
quires a numerical evaluation of the bubble wall solutions 
of the finite T Higgs potential, which is beyond the scope 
of this project. Here, we treat tan/3(T) = V2{T)/v\(T) 
as a free parameter 4 , and we work in the j3(T) <C 1 
limit. Intuitively, we expect j3(T) to be suppressed in 
the limit m 2 H ^> T 2 , since the vev will be confined along 
the {H%) = valley. 

The charge transport dynamics of EWBG are gov- 
erned by a system of Boltzmann equations of the form 
h a = S-f + D a V 2 n a + J2b r °-bn b (39|. Here n a is the 
charge density for species a. The CP-violating source 
generates non-zero n a within the expanding bubble wall, 
at the boundary between broken and unbroken phases, 
due to the spacetime-varying vevs Vi^{T). The diffusion 
constant D a describes how n a is transported ahead of the 
wall into the unbroken phase, where weak sphalerons are 
active. The remaining terms describe inelastic interac- 
tions that convert n a into charge density of other species 
6, with rate r a b. Our setup of the Boltzmann equations 
follows standard methods, described in detail in Ref. (ioj . 

Following Ref. J3j|, we assume a planar bubble wall 
geometry, with velocity v w <C 1 and coordinate z nor- 
mal to the wall. The z > (z < 0) region corresponds 
to the (un)broken phase. We look for steady state solu- 
tions in the rest frame of the wall that only depend on 
z. Therefore, we replace h a — > v w n' a and V 2 n a — > n", 
where prime denotes d/dz. We adopt kink bubble wall 
profiles 

v(T)/T = f [1 + tanh(«/L„)]/(2V2) , (21) 
(3(T) = A/3[l + t&nh(z/L w )}/2 , (22) 

where v(T) 2 = vi(T) 2 + v 2 (T) 2 . We take £ = 1.5, wall 
width L w = 5/T, and T = 100 GeV. Ref. ^ found vi- 
able first-order phase transitions with 1 < £ < 2.5 and 
2 < L W T < 15, depending on the Higgs parameters. For 
dcfiniteness, we take mH 2 = 400 GeV; however, our anal- 
ysis does not account for the crucially important mjj 2 - 
dependence of the bubble profiles. 

Specializing to our 2HDM, the complete set of Boltz- 
mann equations is 

v w n' qa = D q n'g a + 5za{Sf' + TyQ y + T m Q m ) - 2T SS Q SS 

v wn' Ua = D q n u a ~ S 3a(Sf + FyQy + ^mQm) + ^ssQss 

v w n' da = D q n'^ a + T SS Q SS (23) 
v w n' H = Dhu'h + T y Qy - T h Q h 



4 Although the usual tan /3 is not physical at T = 0, the angle /9(T) 
between the T = and T/0 vev directions is physical. 
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with linear combinations of charge densities 




The relevant densities are the ath generation lcft(right)- 
handed quark charges n qa (n Ua , rid a ), and the Higgs 
charge density n B = n Bl + uh 2 (we treat H 12 as mass 
eigenstates in the unbroken phase). We assume that 
(Cabibbo unsuppressed) gauge interactions are in equi- 
librium, as are Higgs interactions that chemically equi- 
librate H12 (provided by Aa^s quartic couplings in V). 
Lepton densities do not get sourced and can be neglected. 
The /c-factors are defined by n a = T 2 k a fi a /6, with chem- 
ical potential [la- 
in the Eqs. (f2"3")l . we take these transport coefficients 
as input: 

Sf w 0.1 x N c \y t y tt \smAtv(T) 2 v w /3(T)'T (26) 
T m a 0.1 x N c \y t vx(T) + y tt v 2 (T)\ 2 T' 1 (27) 

r, » ^a s yfT + 9\y tt fT (^f 2 e^' T (28) 
T ss « Ua A s T , D a « 6/T , D H « 100/T . (29) 

We compute the CP-violating source Sf^ and relaxation 
rate T m , arisin g fo r t^R only, following the vev-insertion 
formalism (4lll42l| (explicit formulae can be found in 
[43j). 5 The sole source of CP violation here is the 
phase 9u = arg(j/ tt ), which is not the same phase that 
enters into i?° s --Bj mixing. 6 The dimensionless nu- 
merical factors (0.1), obtained following Ref. [42l] . arise 
from integrals over tL,R quasi-particle momenta, tak- 
ing as input are the thermal masses (tabulated in (45j ) 
and thermal widths (7^ R a 0.l5g 2 T [H]). The top 
Yukawa rate Ty comes from processes H\ti, O tRg and 
H 2 <->• t R i L HE S3- Tnc strong sphaleron rate T, s [H| 
plays a crucial role in EWBG in the 2HDM [If, dis- 
cussed below, and D q R are the quark and Higgs diffu- 
sion constants [Ho| • The relaxation rate Th is due to Higgs 
charge non-conservation when the vev is non-zero. For 
simplicity, we set = T m [3^]; we find deviations from 
this estimate lead to < 0(1) variations in our computed 
Yb. We have omitted from Eq. (f2"3"|) additional Yukawa 
interactions induced by yt c (e.g., H 2 O tRCr,) because 
we find they have negligible impact on Yr. Moreover, 
CP-violating sources from y tc do not arise at leading or- 
der in vev-insertions. Therefore, yt c plays no role in our 



5 Although there exist more sophisticated treatments, the reliabil- 
ity of quantitative EWBG computations remains an open ques- 
tion (see discussion in [44l0 . 

6 The reparametrization invariant phase is Utt = arg(5/ttj/*i;]'U2), 
but we have adopted a convention where 1)1,2 {T) and yt are real 
and positive. 




-80 -60 -40 -20 



z/K 

FIG. 3: Charge densities in the unbroken phase (z < 0) for 
\ytt\ = 1, 9u = 0.18, v w = 0.05, A/3 = 10~ 2 , giving Y B w 
9 x 10" 11 . 

EWBG setup (this conclusion may not hold beyond the 
vev-insertion formalism) . 

Thus far, we have neglected baryon number violation; 
this is reasonable since the weak sphaleron rate T ws » 
120a^,T [Fl[ is slow and out of equilibrium. Therefore, 
we solve for the total left-handed charge nr, = Ylia n io. 
from Eqs. (|23[) . neglecting T ws , and then treat ni as a 
source for baryon density tir, according to 

v w n' B - D q n B = ~(3T ws n L + TZn B )h , (30) 

with the relaxation rate K = (15/4)r ws [H]. The 
sphaleron profile h{z) governs how T ws turns off in the 
broken phase [53| |. Since the energy of the T = 
sphaleron is i? sp h ~ 4M^ /a w , we take [U 

h(z) = exp(-E sph (T)/T) , E sph {T) = E sph v(T)/v . 

(31) 

Effectively, this cuts off the weak sphaleron rate for rel- 
atively small values of the vev: v(T, z)/T > g 2 /(8ir). 

In Fig. [3l we show the spatial charge densities result- 
ing from a numerical solution to Eqs. (|23p for an example 
choice of parameters giving Y B ss 9 x 10 -11 . In general, 
the individual charge densities have long diffusion tails 
into the unbroken phase (z<0). However, til is strongly 
localized near the bubble wall (z = 0), due to strong 
sphalerons, thereby suppressing tir [49j . This effect can 
be understood as follows: at the level of Eqs. (|2"3"|) . B 
is conserved, implying J2 a ( n qa + "«« + n <i a ) = 0; addi- 
tionally, strong sphalerons relax the linear combination 
of densities 

Q ss w (1/JVc) X]( n ?<. ~ n «o ~ n d a ) ( 32 ) 

a 

to zero. These considerations imply that nr, ~ if strong 
sphalerons are in equilibrium. In Fig. [3J we see that 
strong sphalerons are equilibrated and vanishes for 
z < — 10L W . Since til is non-zero only near the wall, it is 
important to treat the weak sphaleron profile accurately 
in this region, rather than with a simple step function. 
Nevertheless, despite this suppression, EWBG can ac- 
count for Y B hs . (We also note the significant Higgs charge 



7 




0.001 0.002 0.005 0.010 0.020 

v w (wall velocity) 
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FIG. 



4: Computed baryon asymmetry Yb, normalized with 
9 x 10" 11 and (sin M A/3), as a function of 
v w and \ytt\- Vertical axis shows (sinS tt A/3) _1 required for 
viable EWBG. 



respect to Y£ bs 



iih in the broken phase. Although we neglect lepton 
Yukawas here, it is possible that tih could be efficiently 
transferee! into left-handed lepton charge via y^ , thereby 
driving EWBG without suffering from strong sphaleron 
suppression, analogous to Ref. [45|.) 

In Fig. 21 we show how large Yb can be in our model. 
The most important parameters are A/3, y tt , and v w (we 
find Yb is not strongly sensitive to L w or £) . The vertical 
axis shows the (inverse) value of A/3 x sin(9 t t required for 
successful EWBG (Y B = Y° hs ), for different values of 
\yu\ and v w . Our main conclusion is that our model can 
easily account for the baryon asymmetry of the universe 
- even if A/3 is as small as 10 —10 , provided the NP 
Yukawa coupling has magnitude \y tt \ > 0.2, with 0(1) 



phase. Moreover, \y tt \ 
with flavor observables. 



1 is prefered by consistency 



V. CONCLUSIONS 

The dimuon asymmetry reported by DO and the 
branching ratio BR(£> — > tv) [J 0] seem to disfavor the 
CKM paradigm of CP violation in the SM at the ~ 3a 
level. Although more experimental scrutiny is required, 
taken at face value, these anomalies can be accounted 
for by new physics in both B®-B® and B^-B® mixing Q. 
Such new physics would involve new weak-scale bosonic 
degrees of freedom and new large CP-violating phases. 
These two ingredients are precisely what is required for 
viable electroweak baryogenesis in extensions of the SM. 



We proposed a simple 2HDM that can account for 
these B meson anomalies and the baryon asymmetry. 
An interesting feature of our setup is a top-charm flavor- 
violating Yukawa coupling of the new physics Higgs dou- 
blet. The large relative phase of this coupling can explain 
both the dimuon asymmetry and tension in BR(i3 — > tv) . 
Although top-charm flavor violation can give potentially 
large contributions to b —> sj and ek (i.e., less CKM- 
supprcssed than SM contributions), these bounds are 
weakened in precisely the same region of parameter space 



consistent with B^ ■ 



B® s observables. 



We also discussed electroweak baryogenesis. We 
showed that, provided a strong first-order eletroweak 
phase transition occurs, our model can easily explain the 
observed baryon asymmetry of the Universe. CP viola- 
tion during the phase transition is provided by the rela- 
tive phase in the flavor-diagonal th-^R Yukawa coupling 
yu to the new Higgs, and the relevant phase is not re- 
lated to the top-charm CP phase entering flavor observ- 
ables. However, flavor observables and baryogenesis both 
require \y tt \ ~ 1. Additionally, baryon generation is de- 
pendent on a parameter A/3 related to the shift in the 
ratio of Higgs vevs across the bubble wall. We expect 
A/3 to be suppressed in the limit m H ± ^> mw- However, 
we showed that the charged Higgs state can be light 
(m H ± ~ 100 GeV) without conflicting with flavor ob- 
servables due to the large top-charm phase in our model 
(as opposed to the limit m#± > 315 GeV from b — s- sj 
in a type-II 2HDM [Hill). 

It would be interesting to explore the consequences of 
our model for Higgs- and top-related CP-violating and 
flavor-violating observables measurable in colliders, and 
also for rare decays such as K — > irisD. Additionally, 
a more robust analysis of EWBG requires an analysis 
of the finite temperature effective potential in a Type- 
Ill 2HDM, addressing the phase transition strength and 
bubble wall profiles. 
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